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O \ 

t We consider the possibility of formation of an unconventional spin density wave (USDW) in quasi- 

one dimensional electronic systems. In analogy with unconventional superconductivity, we develop 
a mean field theory of SDW allowing for the momentum dependent gap A(k) on the Fermi surface. 
Conditions for the appearence of such a low temperature phase are investigated. The excitation 
spectrum and basic thermodynamic properties of the model are found to be very similar to those of 
d-wave superconductors in spite of the different topology of their Fermi surfaces. Several correlation 
functions are calculated, and the frequency dependent conductivity is evaluated for various gap 
functions. The latter is found to reflect the maximum gap value, however with no sharp onset for 
absorbtion. 
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I. INTRODUCTION 



As a result of intense research during the past few decades, much is known about the properties of density waves 
possessing a constant A order parameter Q. On the other hand it is well known that wavevector dependent order 
' parameters A(k) taking different values at different points on the Fermi surface play an important role in theories of 
superconductivity (termed unconventional) j2| in various systems including the high-T c cuprates j3j. Therefore the 
need to work out the theory of unconventional density waves follows naturally from these precedents. 

This project however is not just of academic interest. The anomalously small magnetic moment in URri2Si2 [Q was 
■ suggested to be explained among other possibilities [|||] by unconventional SDW on a quasi-two dimensional (square) 
Fermi surface JjJ . The basis of this suggestion is the fact that an USDW is in fact not a spin density wave at all, since 
t-H , it is not accompanied by periodic spin density modulation. The order parameter of the phase transition is not the 
spin density itself, but another well defined quantity related to the spin density much the same way as the "effective 
density" responsible for electronic Raman scattering is related to the density operator As another signature 

of an unconventional density wave (DW), low energy excitations due to a DW gap vanishing on certain subsets of 
the Fermi surface may be responsible for the absence of a clear optical gap in the reflectivity data of some of the 
heavy-electron materials |Io|| , including URu2Si2- Optical data with similar features in the low temperature phase of 
quasi-one dimensional M04O11 indicate again the possibility of an unconventional charge density wave (UCDW) state 

r- 00. 

Other candidates for unconventional DW states include the organic conductor a-(BEDT-TTF)2KHg(SCN)4, where 
in spite of a clear indication for a phase transition in magnetotransport measurements neither charge nor magnetic 
order has been established [G~3 T5|, the transition-metal dichalcogenide 2iJ-TaSe2 with its momentum dependent 



CDW gap inferred from angle resolved photoemission studies ||lq|, and the tetrachalcogenide (TaSe^I, for which 



the magnetic susceptibility above the conventional CDW transition temperature shows pseudogap behavior [jl? 18 
without any observable long range charge order. Indeed, a recent suggestion to understand the pseudogap phase of 
the cuprate superconductors also invokes the existence of an unconventional density wave (UDW) of d-symmetry @. 



X 
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Motivated partly by the rich phase diagram of the high T c cuprates, significant steps have already been made towards 
the theory of UDW in a two dimensional electron system |2(J , typically on a square lattice (2^] . 

The objective of the present paper is to develop a detailed theory of UDW-s in quasi-one dimensional interacting 
electron systems. Clearly, the topology of the Fermi surface is radically different in this problem than in previous 
treatments, moreover the strong anisotropy of transport properties in, and perpendicular to the linear chain direction 
is of particular interest. A preliminary report of some of our results has already been published p2fl . The article is 
organized as follows: in Sec. II. we define our model, develop its thermodynamics in mean field theory, and determine 
conditions for the appearance of USDW. Sec. III. is devoted to the calculation of the most important correlation 
functions of the system, while in Sec. IV. we pay particular attention to the optical conductivity of the model. Our 
conclusions are given in Sec. V. 
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II. THERMODYNAMICS OF THE MODEL 



We start with a quasi-one dimensional interacting electron system described by the following one band Hamiltonian: 

H = ^2 £ ( k ) a t* a ^ + ^ H ^( k ' k ^)4+ q , CT a k, CT a+_ qCT ,a kV ,, (1) 
kff k,k',q 
a, a 1 

where aJ CT and ak,cr are creation and annihilation operators of an electron of momentum k and spin c, V is the volume 
of the sample and the kinetic energy spectrum on an orthorombic lattice 

e(k) = — 2t a cos(k x a) — 2% cos(k y b) — 2t c cos(k z c) — [i (2) 

is highly anisotropic (t a 3> tb,t c ). The interaction potential matrix elements in (|l|) are: 

iy(k, k', q) = J d 3 r J d 3 rVk +q (r)^- q (r')^(r - r'V k < (i>k(r), (3) 

where (p^ is Bloch function. In Wannier basis <^k(r) = -^7= e lkR ^(r — R), where N is the number of cells and 
tp(r) is the corresponding Wannier function assumed to be real and an eigenfunction of parity. We note here that 
in a tight-binding solid the Wannier function is well localized, leading to a significant dependence of the interaction 
matrix element (|3|) on the incoming electron momenta k and k'. This turns out to be crucial in order to form an 
UDW, and is readily seen from the expansion including on site and nearest neighbor two center integrals: 

— V(k, k', q, <7, a')= S^ a ^{U + ^[2^ cosgA + 2 J, cos(^ - k[ + qi )Si + 2Fi cos(fc; + h)Si + 

i 

+2C I (cos kiSi + cos k[8 l + cos(fc- - q^Si + cos(fcj + q^S,)]} + 

+5 a: a' y^jVi - Ji)(cOSq l S l - COs(fc; - fc ■ + q l )S i ). (4) 

i 

In the above formula the antisymmetrized (therefore spin dependent) interaction is given with i = x,y,z and Si = 
a, 6, c. The (at most) two center integrals are the on site Hubbard repulsion and the nearest neighbor direct, exchange, 
pair-hopping and bond-charge terms 

U = J d 3 r J rf 3 r>(r)|V(r - r')|^(r')| 2 , (5a) 

Vi = J d 3 r J d 3 r'\^(r)\ 2 V(r - r')\<p(r' - e,)| 2 , (5b) 

Ji = I d 3 r J d 3 rV(r)(^(r' - e t )V(r - r>(r>(r - e*), (5c) 

F t = J d 3 r J d*r'(p{r)<p{r') V (r - r>(r' - e t )^(r - 6i ), (5d) 

d = fd 3 r f d 3 r^(rMr')V(r - r>(r>(r - e*), (5e) 



where ej is the lattice vector in the i direction. 

Although due to its rich structure the interaction in Eq. (|4|) is able to support a variety of low temperature phases 
p3| depending on the Hubbard integrals in Eq.(|^), we are now interested in constructing the mean field theory of an 
USDW. The best nesting vector for the spectrum (0) is obviously Q = (2kp , n/b, tt/c) with the Fermi wavenumber 
kp satisfying /i = — 2i a cos(afcF). In a DW expectation values of the type < a k ~ cr cik+Q,er > will not vanish, defining 

the order parameter of the low temperature phase A(k, a) = | A(k, cr)|e l< ^ k, ' T ) as 

i " 

A(k, a) = - ^( k '> k > Q' o>) < a+ a , a k , +Q ^ > . (6) 

k' <j' 
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Here the overline indicates complex conjugation and the double prime over the summation sign restricts k x values to 
the interval from —2kp to 0. Then the mean field Hamiltonian is diagonalized in the usual way, giving rise to a two 
band quasiparticle spectrum over the new Brillouin zone (0 < k x < 2kp) given by 



g±(k) = g (k) +e (k-Q) ± /e(k)-e(k-Q) (?) 



2 

The new (effectively noninteracting) fermionic quasiparticles are expressed by the original electrons as 

d+ko- = e - *^ k,< ^tt(k, a)a^ a + v(k, cr)a k -Q, -, 

d-^a = e-^'^vik, (7)a kjCT - u(k, <7)a k _ Q>CT) (8) 

with 



u(k, cr) _ 
v(k,a) 




e(k)-e(k-Q) 
2 



Y v v( £(k) 'f' Q) ) 2 +iA(k,cT)p, 

Eqs.(^|) and (0) lead to a self consistency condition for the order parameter known as the gap equation 

k,a V 2 ^ ( s(k)-e(k-Q) )2 + )A(k; 



(9) 



(10) 



where the prime indicates that the k sum runs over the new Brillouin zone only, and f(E) is the Fermi function. 

In the followings we suppress the spin index of the order parameter, since in order to describe an SDW with 
polarization vector parallel to the z axis of the spin space, we can utilize the relation A(k, +) = — A(k, — ) = A(k). 
Moreover, the structure of the gap equation makes it clear that the relevant wavenumber values in the arguments of 
both the gap and the interaction are confined to a narrow region near the Fermi sheet at +kp. Therefore the gap 
equation is in fact an integral equation on the k x = kF plane of variables k y and k z only. Making use of the electron 
spectrum (^) linearized in k x around kp\ £(k) = vp(k x — kp) — 2tt,cos(k y b) — 2t c cos(k z c), we obtain a simplified gap 
equation (vp = 2at a sm(akp) is the Fermi velocity): 







r /c dk z r pkF 


l-TT/b 


2tt . 


I-tt/c 27r JO 



2irv F ^ 2 + |A(k)| 2 

(11) 



Here E(k) = ^[^(k)] 2 + |A(k)| 2 , (3 = l/fc^T, we have neglected terms of order (if,. c /t a ) 2 in the second expression, 
and the kernel of the integral equation is 

P(k WPP P P P 

— ^p— = + cos(kyb) cos(l y b) + — sin(kyb) sm(l y b) + — cos(fc z c) cos(l z c) + — sin(fc z c) sin(^c), (12) 
with coefficients given by 

P = U - V y - V z - J y - J z + (V x + J x ){cos(2k F a) + 1) + 2F X + 8C X cos(fc F a), (13a) 

Pi = -2F y + Jy + Vy, (13b) 

P 2 = 2Fy + Jy+Vy, (13c) 

P 3 = -2P Z + J z + V z , (13d) 

P A = 2F Z + J z + V z . (13e) 



As is seen above in Eq.([12|), the kernel turns out to be diagonal on the basis of the leading harmonics on the {ky, k z ) 
plane. Consequently, the gap will be of the form 

A(k) = A + Ai cos(kyb) + A 2 sin(k y b) + A 3 cos(fc z c) + A 4 sin(fc z c). (14) 
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For vanishing order parameter the five components in ( |11| ) decouple completely, and the critical temperature for the 
development of each type of gap can easily be evaluated. For the conventional SDW with constant gap we recover 
the well known result 

fcB Tj°) = ' 2 LLv F k F e- 2 ' p ^\ (15) 

with 7 = 1.781, the Euler constant, and po(0) = a/irvp is the electron density of states in the normal state per spin 
at the Fermi surface. The critical temperature for the four kinds of unconventional gap formation (j = 1, ...,4) is 

kBT U) = h VFkpe -^po(o)_ (16) 

On cooling the system, that type of an SDW will develop first, for which the T c is the highest. The condition for the 
formation of an USDW of type j is Pj > 2Pq. For example in case of a half filled band an unconventional phase with 
the gap function A(k) = A 2 sm(bk y ) will form if %(V y + J y ) + F y + V Z + J z > U + 2F X . Clearly, a combination of 
interchain Coulomb and exchange integrals overwhelming the on site repulsion will facilitate the development of an 
USDW (negative interchain pair-hopping integral favours a gap with cosine dependence). 
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FIG. 2. The density of states for conventional (full line), and unconventional (dotted line) density waves. 
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The quasiparticle energy spectrum from Eq. (|7|) is shown in Fig. [l| for the two typical unconventional gap functions 
of k y (we neglected any dispersion in k z here for clarity) . The excitation energy vanishes on lines (note the additional 
k z direction in the Brillouin zone) of the Fermi surface, and this will determine the nature of the thermodynamics of 
the system. The corresponding density of states (DOS) is calculated as 

p{e) = r d(bk y ) r d( C k z ) \e\ 

Po(0) y_ w 2tt J_„ 2ir ^-|A(k)|a' 

and is shown in Fig. ^| for both the well known conventional situation and for the unconventional cases determined 
analytically by p(E)/ Po (0) = {2\E\/n\A j \)K(\E\/\A j \) if \E\ < |A,|, and p(E)/p o (0) = (2/tt)K(\A,\/\E\) if \E\ > 
\Aj\. In the latter case the DOS vanishes linearly at the Fermi energy and diverges logarithmically at the maximum 
gap value, as follows from the properties of the complete elliptic function of the first kind. 




T/T c 

FIG. 3. The temperature dependence of the gap amplitude for conventional (Ao) and unconventional (Aj) density wave. 

Assuming that only one kind of gap (cither the conventional or one of the four unconvcntionals, whichever opens 
at the highest T c ) persists all the way down to zero temperature, we can use Eq.(|ll|) in order to evaluate the 
temperature dependence of the gap amplitude. This is the same for all unconventional gap types, and is shown on 
Fig. U along with the conventional dependence, displaying only small differences between the two. At zero temperature 
the unconventional gap takes the value |Aj(0)| = (2n /jy/e)kBTc : '\ leading to a gap maximum to T c ratio of 4.28, 
instead of the ratio 3.52 in the conventional case. For T ^T c the unconventional gap decreases from its T = value 
as |Aj(T)/Aj(0)| = 1 — 3('(3)[A:BT/|A : j(0)|] 3 , which is to be contrasted with the exponential correction for a constant 
gap. Similar power law behavior is found for other quantities as well, due to the line nodes in the excitation spectrum. 
For example the specific heat vanishes like T 2 at low temperatures, and normalized to the normal state value we obtain 
Cv\T)/d n \T) = [27C(3)/7t 2 ][A:bT/|A : ,(0)|], as opposed to the exponential freezing out for a conventional gap. Close 

to the transition temperature the unconventional gap vanishes like [|A J (T)|/fc B T c (: ' ) ] 2 = [32tt 2 /21C(3)][1 -T/T^}, i.e. 
in a square root manner with somewhat different prefactor compared to the conventional case (87r 2 /7C(3)). The mean 
field transition leads to a specific heat jump at T c with the relative value of ACi ' /Cy — 8/7£(3) = 0.95, exactly 
two third of the conventional value. Due to the presence of line nodes in the gap function the above thermodynamic 
properties are identical to those of a d x 2_ y 2 superconductor [^) in spite of the different topology of their Fermi surfaces. 



III. CORRELATION FUNCTIONS AND THE NATURE OF THE ORDER PARAMETER 

At this point it is important to call the attention to the fact that while in a conventional SDW the order parameter 
is proportional to the Fourier component of the magnetization density at the nesting vector Q, an unconventional 
order parameter has nothing to do with the magnetization. In order to see this we evaluate the magnetization using 
the transformation in Eq.(||): 
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< m(Q) >= 



-(IB ^ (T < a 

k,(T 



k-Q,<r 



k 



A(k) tanh{|£:(k)} 
W) ' 



(18) 



It is easily seen that Eq.(|i8|) yields zero magnetization for any of the four unconventional gap functions, as opposed 
to the conventional situation leading to < m(Q) >= 2N/j,b^o/Pq- This means that an USDW is not accompanied 
by a spatially periodic modulation of the spin density, although the expectation value < a^Q ^ak.o- > becomes finite 
and the existence of a robust thermodynamic phase transition is unquestionable. This feature of the unconventional 
density waves makes them suitable candidates for explaining low temperature phase transitions, where conventional 
order parameters such as charge-, or spin-density modulation are not observable, like in a-(BEDT-TTF)2KHg(SCN)4 
p3| , or in URu2Si2 Q respectively. 

What kind of physical quantity does then an unconventional gap correspond to? Again utilizing the gap equation 
(pi]), we can convince ourselves that if for example an unconventional gap Ai cos(k y b) develops at low temperature, 
then the quantity 



&(q) 



asm[b{k y + q y /2)]a£ „a k+q , CT 



(19) 



assumes a finite expectation value < S^Q) >= NA\/Pi. Therefore the Fourier component with wavenumbcr Q of the 
"effective" spin density S z (r) plays the role of the order parameter in this unconventional case. Clearly, experimental 
observation of this order parameter is possible only in probes coupling directly to this physical quantity. The situation 
is rather similar to electronic Raman scattering, where the photon-electron vertex carries momentum dependence, 
and measures the effective density correlation function ||p|, instead of scattering on just density fluctuations. 

The rest of this chapter will be devoted to the evaluation of certain correlation functions which are of particular 
interest. We begin with the charge susceptibility Xnn(<l,t) = i < [n(q, t), n(— q, 0)] >, the autocorrelation function 



of the density operator n(q) = 



k.cr k,(7 



«k+q,CT- The quasiparticle contribution to the frequency dependent charge 



susceptibility in the long wavelength limit reads as 



, , | e(k)g(k + q)+Re(A(k)A(k + q)) ^ rj , r ^ u 
X „„(q, W )= -^^1+ E (k)E(k + q ) ) + q 



)]-f[E(k)]}x 



lu + E{k) - E(k + q) lu - E(k) + E{k + q 

( 




£(k)g(k + q)+Rc(A(k)A(k + q)) 
E(k)E(k + q) 



x{l-f[E(k + q)]-f[E(k)}} 



\ uj + E(k) + E(k + q) uj-E(k)-E(k + q) 
while for wavenumbers around the nesting vector we obtain 



(20) 



Xnn(Q + q,w) 



k.tr 



f[E(k + q)]-f[E(k)] ( 
uu + E(k) - E(k + q) \ 

f[E(k)]-f[E(k + q)] 
uu - E{k) + E(k + q) 

l-f[E(k)]-f[E(k + q)] 
v + E(k) + E(k + q) 

f[E(k)] + f[E(k + q)]-l 
lu - E(k) - E(k + q) 




£(k + q) 



E(k^ 


-q)y 




4) 


S(k4 




£(k4 


q)\ 


E(k^ 








E(k 


+ q) / 



+ 



+ 



(21) 



The analytic structure of the spin susceptibilities are very similar to Eqs.(|20|) and (|2lJ). For example the longitudinal 
spin susceptibility xs~s z — Xnn/^ for all wavelengths, and the transverse spin susceptibility Xs+s- — Xnn/^ for short 
wavelengths. For long wavelengths however, we encounter coherence factors different from those in Eq.(pQ), namely 
the sign of the two Re(A(k) A(k + q)) terms becomes negative. 

In the followings we evaluate the above mentioned correlation functions in both the static (first lu — > 0, then q — > 0), 
and the dynamic (first q — > 0, then lu — > 0) limit. 



G 



A. Susceptibilities in the long wavelength limit 
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FIG. 4. The temperature dependence of the long wavelength static charge susceptibility Xnn f° r conventional (dotted line) 
and for unconventional (dashed line) density waves. 
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FIG. 5. The temperature dependence of the long wavelength dynamic transverse spin susceptibility Xg+g- f° r conventional 



(dotted line) and for unconventional (dashed line) density waves. 



The dynamic limit of the long wavelength charge susceptibility is trivially zero, therefore we begin with the static 
limit of Eq.©; 



\ 0S 

Ann 



V 



1 I'OC 

J2{-f'[E(k)}} = ffo(O) / dE[-f(E)}p(E)/p o (0), 

I. _ J —oo 



(22) 



where f'(E) is the derivative of the Fermi function, and <?n(0) = 2iVpo(0)/V" is the normal state DOS (per unit 
volume) at the Fermi energy including spin degeneracy. The above susceptibility takes this value for T > T c , and its 
temperature dependence is shown in Fig. || for both the conventional and the unconventional cases. We note here, that 
following superconductivity terminology, Eq.(|22|) is often expressed as \nn = ffo(0)(l — fs)> where / s is the superfluid 
fraction. In our case of course, condensate fraction is more appropriate. Due to the above mentioned relations between 
correlation functions, Fig. |] represents also the longitudinal spin susceptibility, experimentally accessible through 
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Knight-shift measurement, vanishing at low temperatures. This is in contrast to the static homogeneous transverse 
spin susceptibility X^+s- = ^ X J2k I A(k)| 2 /[i?(k)] 3 = ffp(0)/2, which is independent of temperature. Returning 
to the charge (or longitudinal spin) susceptibility, Eq.(|2|) is evaluated at low temperatures yielding exponential 
freezeout for conventional SDW, and x°n/ffo(0) = 2 m(2)/c B T/| Aj(0)| for the USDW. Close to T c on the other hand 

xl S n/9o(0) = (4T/TP - l)/3 for USDW (instead of the conventional value 2T/T C (0) - 1). 

Wrapping up our discussion of the long wavelength correlation functions, we consider the dynamic limit of the 
transverse spin susceptibility 



,0D 



Xs+s- 



- £ i^Cu - mm}} - ^ f dE l ~ 2f } E) r « f dick ^~ - |A(k)l (23) 



V^[E(k)}^ JL v JU 2 Jo ' £ 2 J_„ 2tt J_ n 2vr - |A(k 



which is nonzero in the DW state, as shown in Fig. ^ for both conventional and unconventional situations. The 
susceptibility normalized to its zero temperature value Xs+s-(T , )/[ffo(0)/2] = fd is again the condensate fraction, 
but now in the dynamic limit. In order to see the difference compared to the static condensate fraction, we first 
realize that [g (0)/A}{f d - f s ) = £jj-f'[£;(k)]}| A(k)| 2 /[£(k)] 2 , then consider the limiting values for low and 
high temperatures. For T — > the dynamic condensate fraction /<j = 1 — ln(2)fcsT/|Aj(0)| in USDW (the finite 
temperature correction is exponentially small for the conventional case). Close to T c the dynamic condensate fraction 
vanishes like 7r|Ao|/4fcaT for constant gap, and like |Aj|/2fceT for momentum dependent gap. 



B. Susceptibilities around the nesting vector 



We turn our attention now to the behavior of the short wavelength correlation functions based on Eq.(|2l|). These 
are of particular interest in describing phenomena related to the phase transition involving the nesting vector Q. It 
is easily seen, that in the normal state (T > T c ) the charge susceptibility at the nesting vector follows a logarithmic 
temperature dependence 



-t 



i - 2/K(k)] 

2£(k) 



0o (0) . 2"fv F k F 

In 

2 irk B T 



(24) 



for both conventional and unconventional cases independent of the limiting procedure (static or dynamic). 
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FIG. 6. The temperature dependence of the short wavelength charge susceptibility \nn in the static limit for conventional 
(dotted line) and for unconventional (dashed line) density waves. We have chosen PoPo(O) = Pjpo(0)/2 ps 0.43. 

Below the transition temperature we obtain a still rather simple formula 



X Q 

Ann 



1 A l-2/[i?(k)] g (0) / 



V 



2E(k) 



(25) 
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incorporating all the variations in the condensate fraction /, discussed extensively in the previous subsection in 
all limits and cases. It is easy to deal with the first term in Eq.(^H|) for a constant gap, since the gap equation 
(O) relates just such an expression with the inverse of the coupling constant. Therefore in the conventional case 
X% n = 2N/VPq — [<7o(0)/4]/, leading to a monotonically decreasing susceptibility in both the static and dynamic limit 
as the temperature is lowered (see Figs. || and 0). 
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FIG. 7. The temperature dependence of the short wavelength charge susceptibility \nn in the dynamic limit for conventional 
(dotted line) and for unconventional (dashed line) density waves. Coupling constants are chosen as in Fig. ^| 

In case of a momentum dependent gap, the first term in Eq. (p5|) is not exactly what appears in the gap equation 
(|ll|), but at the expense of a correction factor we can bring in the inverse coupling constant as in the conventional 
case: X % n = ^N/VP, - [ffo(0)/4](/ - Sf), where 



Sf = 4 



d{bk y ) r d{ck z ) 



2tt 



2tt 



dERe- 



1 - 2f(E) 

y/E*-\A(k)\< 



A(k) 



(26) 



This correction factor is evaluated as Sf = 1 — 41n(2)fcsT/|A. ) (0)| for low temperatures, and as Sf = (2/3)(l— T/T c ) 
close to the critical temperature. 

Our results concerning the susceptibility at the nesting vector are summarized in Figs. ^ and [t] for the static 
and dynamic limits respectively. In the conventional situation Xnn °^ Xs s 1S P e& ked at the critical temperature 
signaling the phase transition. Indeed, if we consider the full spin susceptibility in random phase approximation 
(RPA), the Stoner denominator vanishes exactly if we approach from above, leading to divergent response. The 
other (unconventional) coupling channels do not contribute to the charge (spin) response due to their momentum 
dependence. In the unconventional case however, RPA corrections will not lead to divergence, since the dominant 
unconventional channel does not couple to charge or spin density, while the conventional coupling constant is not 
strong enough to make the Stoner denominator vanish. Instead, the autocorrelation function of the effective spin 
density S z (see eg. Eq.(|l9])) will be divergent at Tc in RPA, as it should if we are to have an unconventional phase 
transition. 



IV. OPTICAL CONDUCTIVITY 



The frequency dependent conductivity provides a wealth of information about both the quasiparticle and the 
collective excitation spectrum of density wave materials |IJ] . Therefore in this section we investigate the properties of 
the conductivity tensor a a p{uj) = K a p(q — Q,Lo)Jiui in our model of USDW. The electromagnetic kernel consists of so 
called paramagnetic and diamagnetic parts: K a p = K^p + K^p- The diamagnetic term = — e 2 V~ x a m^(k) 
involves the effective mass tensor derived from the electronic spectrum (^|). The paramagnetic part K^g = Xj a jp 1S the- 
correlation function of the corresponding components of the current operator given by j (q) = — e v(k)a^ ^ak+q,^ 
in the q — * limit, where the electron velocity is again obtained from Eq.(0). The current correlation function turns 
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out to be given by the same formula as in Eq. (pp|) , except for the different coherence factors and the multiplicative 
term e 2 i> Q (k)u^(k) under the summation. Since our system is ideal in a sense that it does not include any source of 
damping for electrons (for example impurity scattering), the real part of the conductivity consists of a sharp Drude 
peak and a regular contribution: Re<7 a| g(w) = irD a pd(uj) + Recr^ 9 (w), where the Drude weight is given by 

(27) 



D °e = V £ { m ^ (k) ~ Va{k)Mk)l ^W {1 2f[m]} 



The Drude peak is the only component in the normal state, but its weight decreases below T c and vanishes completely 
at zero temperature. For example in the chain direction the Drude wight is related to the dynamic condensate 
fraction by D xx — e go(0)v F (l — fd). The missing oscillator strength is taken over by the regular component at finite 
frequencies: 

Rea-M = 9o( 0)4tanh ( ^L) f « f ^ Re ^(k^g 

af)K ' y K u 2 \4k B Tj J_„ 2tt J_ n 2tt ^( w /2) 2 - |A(k)| 2 V ' 

In order to obtain characteristic lineshapes for the optical conductivity tensor, we first realize that it is diagonal, 
then consider various gap functions and electric field directions. In case of a conventional SDW the conductivity is 
given by 

B*r> > 0) . «W. tani ) !l!M Re _i_, (29) 

where v x = Vf, v y — \[2bt\ ) and v z — \/2ct c . The frequency dependence of the conductivity is the same for all three 
directions of the electric field, and is shown in Fig. || at zero temperature. We recognize the sharp onset of absorbtion 
at lu = 2|A | due to the constant gap. 
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w/|Ao| 

FIG. 8. The real part of the complex conductivity of a conventional SDW at T — 0. 

In the unconventional situations we expect finite absorbtion below the maximum gap, since we have low energy 
optical excitations around the nodes of the order parameter. Consider first the conductivity in the chain direction 

Re<r (" > 0) = e 2 5o(0)4 tanh ^M/ , (30) 

where I(k) = f* /2 dtp sin 2 ^Re(l - fc 2 sin 2 (p)~ 1/2 , evaluated as I(k < 1) = [K(k) - E(k)]/k 2 , and I(k > 1) = 
[K(X/k) — E(\/k)]/k. This function is plotted in Fig. ^| at zero temperature. We see a logarithmic divergence 
at the maximum optical gap, and a substantial absorbtion below that gap as expected. We note here, that if we 
consider the conductivity in one of the directions perpendicular to the chain, we obtain the same lineshape if the 
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unconventional gap varies in the other perpendicular direction. Of course we need to replace vf in Eq.(|30|) by the 
proper perpendicular velocity. At this point it is appropriate to remark that we are calculating only quasiparticle 
contributions to the conductivity, although it is known that due to the sliding motion of the condensate a collective 
contribution will also be observed in the chain direction [Q. However, if the electric field is perpendicular to the 
chains, no such contribution is expected, therefore our results are directly applicable to the experimental situation. 



b 




FIG. 9. The real part of the complex conductivity of an USDW in the chain direction at T = 0. 




FIG. 10. The real part of the complex conductivity of an USDW in the y direction at T = 0, A(k) = A2 sm(k y b). 



Finally, we consider the situation when the electric field is aligned in that perpendicular direction in which the 
order parameter varies. Without loss of generality we can take this to be the y direction. There are two possibilities 
corresponding to the gap function A(k) = AiCOs(fc y 6), and A(k) = A2 sin(/c y 6). Eq.(p8|) leads to the following 
expressions for the conductivities: 



and 



Rca c ° s (u; > 0) = e 2 5o (0)^tanh 



Reo-,%"(w > 0) = e 2 5o (0)^tanh 



4fc B T 



4fc B T 



|Ai| 



Ic, 



8|A 2 | 2 



2 * sin 



2|Ai| 



2IA, 



(31) 



(32) 
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where I S in{k) — J d(/?sin 4 y>Re(l — k 2 sin 2 ip) 1 / 2 , and I cos {k) — I(k) — I S i n (k). The former function is evaluated 
as I sin (k < 1) = [(2 + k 2 )K{h) - 2(1 + k 2 )E(k)]/3k 4 , and I sm (k > 1) = [(1 + 2k 2 )K(l/k) - 2(1 + fc 2 )£(l/fc)]/3fc 3 . 
We plot these results in Fig. |l(J and in Fig. [ll] for the sine and cosine dependence of the order parameter respectively. 
As seen on Fig. |l0[ due to the matching of the k dependences of the gap and the velocity v y (]i) = 2btb sm(k y b) the 
low frequency conductivity is suppressed and proportional to u> 2 , while the logarithmic divergence is still there at 
u) = 2 1 A.2 1 , as for the chain direction conductivity (see Fig. ^). This reasoning is quite similar to the one used in 
explaining the lineshape of the B\ g Raman response in d x 2_ y 2 superconductors |25[ . 
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FIG. 11. The real part of the complex conductivity of an USDW in the y direction at T = 0, A(k) = Ai cos(k y b). 

In case of a cosine gap on the other hand, the velocity is zero at the gap maximum, therefore the logarithmic 
divergence is cut off at ui = 2|Ai|, and we have a monotonically decreasing conductivity as shown in Fig. |llj. 

V. CONCLUSIONS 

In this paper we have developed the mean field theory of unconventional density waves in quasi-one dimensional 
interacting electron systems. Although our calculations refer explicitly to unconventional spin density waves (USDW), 
most of the results apply to unconventional charge density waves as well. We have found that the excitation spectrum 
and thermodynamics of our model are identical to those of a d-wave superconductor, due to the presence of line nodes 
in the gap function A(k y ,k z ). Formation of an USDW is facilitated by a combination of interchain Coulomb and 
exchange integrals overwhelming the on site repulsion. 

It is important to realize that an USDW is not accompanied by a spatially periodic modulation of the spin density. 
Instead, an effective spin density plays the role of the order parameter, which is observable only in probes coupling to 
this quantity through form factors significantly dependent on wavenumber. This feature of the unconventional density 
waves makes them suitable candidates for explaining low temperature phase transitions, where conventional order 
parameters such as charge-, or spin-density modulation are not observable. This state of affairs is sometimes referred 
to as " hidden order" . The above conclusion is also corroborated by our investigation of charge and spin correlation 
functions at the nesting vector, which do not diverge at the critical temperature in random phase approximation, 
unlike the effective spin density correlator, which does. The homogeneous static spin susceptibility shows qualitatively 
the same anisotropic behavior below the transition temperature as in the conventional case. 

Finally, we have calculated the quasiparticle contribution to the frequency dependent conductivity for an USDW 
system without quasiparticle damping (collisionlcss limit). The lineshape always exhibits absorbtion for frequencies 
below the maximum optical gap in the quasiparticle spectrum, but varies significantly depending on the direction of 
the electric field and on the functional form of the gap function. These differences can be exploited in determining 
the nature of the condensate by optical spectroscopic tools. 
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